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Abstract. We investigate the short-distance statistics of the local density of states i/ in long one dimensional disor- 
dered systems, which display Anderson localization. It is shown that the probability distribution function P{v) can 
be recovered from the long-distance wavefunction statistics, if one also uses parameters that are irrelevant from the 
perspective of two-parameter scaling theory. 

PACS. 72.15.Rn Localization effects (Anderson or weak localization) - 05.40.-a Fluctuation phenomena, random 
processes, noise, and Brownian motion - 42.25.Dd Wave propagation in random media - 73.20.Fz Weak or Anderson 
localization 



1 Introduction 



Wave localization in a disordered potential is the most strik- 
ing hallmark of systematic interference by multiple coherent 
scattering I1I2I3I4I5I6I7I8I . Systematic constructive interfer- 
ence in a spatially localized region results in a confinement of 
the wavefunction, which decays exponentially away from the 
localization center (with a decay length koc, the localization 
length), in contrast to the extended waves in constant or spa- 
tially periodic potentials. Localization comes along with large 
fluctuations of the wavefunction, which can be induced by chang- 
ing the disorder configuration. The wavefunction statistics can 
be probed, e.g., globally across a system of finite length Lgys 
by the dimensionless conductance (transmission probability) g, 
and inside a semi-infinite system (Lgys = oo) by the local den- 
sity of states 1/ at a distance Lopcn to the opening. 

Theories of localization often focus on the long-distance 
wavefunction statistics, where a high degree of universality pre- 
vails. For instance, distribution functions are restricted to log- 
normal forms as a consequence of the central limit theorem, 
which leads to two-parameter scaling (TPS) |9|. Consequen- 
tially, for the local density of states, the probability distribution 
function P{i^) is characterized by the mean logarithm c['^^ = 

— (In u) and its variance Cj"^"* = var In u. The TPS observation 
has found many applications IIOII 111211 3I14I15I16I . An even 
enhanced degree of universality arises in the random-phase ap- 
proximation (RPA), where single-parameter scaling (SPS) ap- 
plies 1 17 18 19 20 1, and both parameters further are connected, 

e.g.byC^"^) ^2C|'^)for one-dimensional systems |21|22|23J. 
(It was recognized very early that SPS breaks down for strong 
disorder, see, e.g., Ref. |24|.) 

In this paper we point out a connection of the long-distance 
statistics to the short-distance statistics in the one-dimensional 
Anderson model of localization, probed by the local density of 



states i^. Namely, we find that the distribution function P{iy) for 
short distances reliably can be approximated with the help of 
parameters that are extracted from the long-distance limit, in- 



cluding parameters (besides C^^^ and Cj"'') that are irrelevant, 
from the perspective of TPS, for the long-distance wavefunc- 
tion statistics themselves. 

We start this paper by an analysis of P{i^) and P{g) from 
the perspective of large-deviation statistics |25|, which goes 
beyond the central-limit theorem, and identify quantities C„, 
rt > 3, which are irrelevant for the long-distance wavefunction 
statistics, but will turn out to be useful for the short-distance 
wavefunction statistics. Each quantity defines its own length 
scale by its asymptotic slope c„ = limi_,oo dCn/dL (where 
L = Lsys for g and L = Lopon for v), in analogy to the rela- 
tion between Ci ~ 2L/l\oc and and the localization length 1\qc- 
The length scales obtained from v and g coincide. The constant 
offsets dn = lim^^oo Cn — Lcn are shown to contain infor- 
mation on the reflection phase, which allows to test the RPA. 

Then we discuss that P{v) for short distances iopcn ^ ^loc 
can be reconstructed from the parameters c„ and dn- This is in 
striking contrast to P{g), for which the parameters show a tran- 
sient behavior for small Lsys (where wavefunctions are not yet 
localized), as was pointed out very recently in Ref. \ Yii (see 
also Ref. |26|). Our observations for the short-distance statis- 
tics lead us to conclude that the cumulants C„ are useful char- 
acteristics of localized wavefunctions, even though they are not 
relevant in the long-distance Umit because of TPS. 

Finally, we analytically and numerically investigate the pa- 
rameters Cn and dn in various regimes of the one-dimensional 
Anderson model. 

The paper is organized as follows: In Section |2] we de- 
scribe the general implications of large-deviation statistics for 
the scaling of the distribution functions P{g) and P{v), and 
identify the parameters c„ ad dn in the cumulants C„, n > 3. 



2 



H. Schomerus, M. Titov: Short-distance wavefunction statistics in one-dimensional Anderson localization 



In Section |3] we specialize to the one-dimensional Anderson 
model. In order to motivate subsequent considerations, we first 
illustrate in Section lTTI the length dependence of the cumulants 
C„ by numerical simulations. Then (Section]^} we briefly re- 
view the analytical theory for the asymptotic slopes c„ 1 26 27 1 
and extent it to the case of competition between onsite disorder 
and offsite disorder close to the band center. We also present 
the theory for the asymptotic offsets dn- In Section lX3l we in- 
vestigate the dependence of the parameters in various regimes 
of the Anderson model. Our conclusions are given in Sec.|4] 

In order to facilitate a parallel discussion of the statistics of 
g and v, we use the common notation L = Lgys when consid- 
ering g and L = Lopcn when considering i/. One has to bear in 
mind that in the latter case, Lgys = oo and hence one always 
discusses localized wavefunctions, while in the former case this 
is true only for L = Lgys ^ ^loc- 



This paper is centered around our numerical observation in 
Sec. |3]that Eq. ilhi holds even for short distances to the open- 
ing -Lopcn 'loc, and hence can be used in regions where the 
central-limit theorem does not apply. This makes the parame- 
ters c„ and dn with n > 3 observable in the distribution func- 
tion P(j^), while in the long-distance behavior only ci and C2 
are relevant parameters 1 9 1 . 

Presently, analytical results for the distribution function P{g) 
and P(i^) for short distances are only available in the regime 
of single-parameter scaling. The local density of states obeys 
a strict log-normal distribution for all distances |23 30 1, and 
hence complies with our central observation. Equation ( fTal can- 
not be extended to short distances, even in the regime of single- 
parameter scaling 1.21 1 : for studies outside this regime see, e.g., 
Refs. 1131261 . 



2 Large-deviation statistics 

Large-deviation statistics often is introduced as the third and fi- 
nal step in a progressively refined analysis of the asymptotic be- 
havior of probability distribution functions, where the first step 
is the law of large numbers and the second is the central-limit 
theorem. In localization, the law of large numbers certifies that 

the Lyapunov exponent 7 = c'f^ /2L is self-averaging in the 
limit L ^ 00 f281, with asymptotic value lim^^ooT = ^loc- 
The central-limit theorem delivers a statement about the finite- 
length corrections to this asymptotic value, which are char- 
acterized by Cj""* : The variance var 7 — C^'^ / decreases 
asymptotically as L~^. Presently, we find it useful not address 
the Lyapunov exponents, since these are defined with help of 
the system length L, but to rely on quantities that only involve g 

or V, like c[^\ &^^\ and C^''\ The law of large numbers 
and the central-limit theorem predict a linear growth of these 
quantities with L. The full picture is unfolded in the frame- 
work of large-deviation statistics |25|: All cumulants can in- 
crease linearly with length or distance, 

C(f) = (((- ln.g)")) ~ c(f)L + (L » k,,), (la) 

= (((- InuT)) ^ c^:^L + (L » /,oc), (lb) 

where the coefficients dn are the subleading corrections that 
can be neglected in the asymptotic limit, but will be seen to en- 
code information on the reflection phase that allows to test the 
validity of the RPA. For the conductance g, the linear scaling 
of the cumulants C^^'' with L and the connection of the rfi^' to 
reflection phases also has been found in a constructive theory 
by Roberts 1291 . The parameters Cn can be extracted from the 
averages 

c(^) iO = - lim - \n{g-') = ^ i-c^f ) (2) 

n 

(or equivalently for i') as function of the continuous parameter 
^. Note the exponential dependence of the moments on L due 
to localization, in contrast to the power-law dependence in the 
critical regime around a metal-insulator transition jTJ . 



3 One-dimensional Anderson model 

The previous section |2] put forward some very general argu- 
ments from large-deviation statistics. The relevance of the asymp- 
totically defined parameters c„ (and dn), n > 3 for finite- 
distance wavefunction statistics, and the question whether these 
parameters indeed contain information that is independent from 
what is encoded in the parameters ci and C2, only can be an- 
swered by a direct investigation. In the following we will ana- 
lyze the wavefunction statistics in the one-dimensional Ander- 
son model f3l, given by the Schrodinger equation discretized 
on a chain (lattice constant a = 1) 

ii-iV'i-i + iiV-i+i ^{Vi~ E)iPu (3) 

where the hopping matrix elements ti and the disorder poten- 
tial Vi are random. We assume box distributions with {ti) — t, 
{Vi) = 0, iUtn,) = e + \Dt5in,, and {ViV„,) = 2Dv5hn- 
Without any restriction we can set t = 1, which fixes the en- 
ergy scale in the dispersion relation E(k) = — 2cosfc of the 
clean system {Dt = Dy = 0). The disorder strength will be 
characterized by the perturbative mean-free path 1311 

Zpcrt = (4 - E^)/{Dv + Dt), (4) 

and the balance between onsite and offsite disorder will be 
characterized by the parameter 

5^{Dv~Dt)l{Dv + Dt). (5) 

First we will present the results of numerical simulations to 
illustrate the usefulness of the cumulants C„. Next, in order to 
give a flavor for the mechanism behind the asymptotic linear 
growth Q of the cumulants for the specific case of wavefunc- 
tion localization, we extent the analytical theory of Refs. II26I 
|271 for the asymptotic slopes c„ to the case of competition of 
onsite- and offsite disorder close to the band center E = Q, and 
also present the theory for the offsets dn- The extension equips 
us with a means to violate SPS, which is finally compared to 
other means in order to determine the mutual (in)dependence 
of the parameters c„ and dn- 
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Fig. 1. Cumulants Cn calculated from g as a function of L = Lsys, 
and from u,v2&& function of L = Lopen- The data points are obtained 
by a numerical simulation of the one-dimensional Anderson model 
with onsite disorder, for E = Q, Dv = 1/75, (left panels), E = ~1 
, Dv = 1/100 (middle panels), and E = -1.95, Dv = 0.0006583 
(right panels). This corresponds to weak disorder, with a perturbative 
mean-free path of Ipat ~ 300 in all cases [see Eq. ^4). The lines are 
analytical weak-disorder predictions of the asymptotic linear behavior 
C„ ~ c„L, taken from Ref. |251 for iS = and following the RPA 
for£ = -1 andE = -1.95. 



3.1 Numerical illustration of the cumulants C„ 

Here we illustrate the length dependence of the cumulants C„ 
by the results of numerical computations in ensembles of 10^ — 
10* disorder realizations. The results for different strengths of 
onsite disorder (6 = 1) are presented in Fig.n(Zpcrt = 300), 
Fig. EK/port = 24), and Fig. 13 (Zpcrt = 1.5). Three repre- 
sentative values of energy are chosen: (a) E ^ Q at the band 
center, (h) E = -1 in the SPS region, (c) E = -1.95 close 
to the band edge. [The constant perturbative mean free path 
/pert in any figure has been obtained by adjusting the disor- 
der strength according to Eq. (|4j; for values see the figure cap- 
tions.] The significance of these three regions of energy will be 
discussed in the following Subsection 13.21 Plotted as a func- 
tion of length are the cumulants C„ calculated from g and 
ly, as well as from the 'mesoscopic' local density of states i/, 
which is obtained from ly by averaging over a Fermi wave- 
length Af = 27r/ arccos(-£'/2) (with Xp ^ 4 for E = 0, 
Af = 6 for = -1, and Xp w 28 for E = -1.95). The 
mesoscopic density of states accounts for a limited resolution 
that may be encountered in an experiment. It discards the nodes 
of the wavefunction (whose impact strongly depends on the di- 
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Fig. 2. Same as Fig.Q but for stronger disorder with Ipat ~ 24: 
Dv = 1/6 (left), Dv = 1/8 (middle), Dv = 0.00823 (right). The 
lines in the right panels (c) are the analytic predictions for the given 
disorder strength close to the band edge, taken from Ref. 1271 . 



mensionality of the system) and only captures the smoothly 
varying envelope (which is more robust). 

The cumulants all increase linearly for L ^ l\ac, and one 
may associate a length scale limL^oo 2L/C„ = 2/c„ to each 
of them. The slopes c„ are identical for all three underlying 
objects, and hence for the sets of parameters {cn^ } = {cn^ } = 
{c^''-'} = {c„} coincide. 

As advertised above, the cumulants Cji^ and ci^-* increase 
linearly already for small L < l\oc (moreover, the offsets for D 

are vanishingly small), while the cumulants ci^'* become linear 
only after some transient length, see Ref. 1131 (these cumulants 

also have a finite offset rfi^'). This means that the asymptoti- 



cally defined parameters Cn'' and d\i' can be used to estimate 

the short-distance behavior of C^"^^ and Cn'^. In order to esti- 
mate the distribution functions, parameters with n > 3 have to 
be included, since the central-limit theorem does not yet apply 
for short distances. This is displayed in Fig.|2 which compares 
P(lnj/) for L = l\oc/% E — —1, Zport — 24 with a normal 



distribution, which only accounts for (7}"^^ and 02"', and with 
a generaUzed normal distribution (the so-called Pearson system 



P{x) = C(a + &.T + cx2)-i/2^ 



X exp 



(6 -I- 2cm) arctan[(6 + 2cx) /\/4:ac - 6^] 



62 



(6) 
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Fig. 3. Same as Figs.Qand|2| but for stronger disorder with Zpert = 
1.5: Dv = 8/3 (left), Dy = 2 (middle), Dy = 0.1316 (right). The 
lines in the right panels (c) are the analytic predictions for the given 
disorder strength close to the band edge, taken from Ref. 1271 . 



which accounts for the first four cumulants by the four con- 
stants a, b, c, and to. The cumulants have been reconstructed 
from their asymptotics ilh\ (hence, from the asymptotically de- 
fined quantities On ^ and rfi'^^ ), and differ from the numerical 
values of the data by less than three percent. 



3.2 Analytical theory for the slopes c„ and offsets dn 



3.2.1 Slopes c„ 



Recently I26I27I . we have been able to extent Halperin's phase 
formalism I33I22I . which allows to calculate l\oc and hence 
ci, to all slopes c„. This formalism can be applied for arbi- 
trary Xp/lpcrt, i-e., also for relatively strong disorder, as long as 
^pert 3> 1 (the lattice constant, set to unity in this paper). Other 
formalisms like the supersymmetric a model and the Berezin- 
skii technique are rather more restrictive and cannot directly 
address the logarithm of g and i'. It turned out that three dif- 
ferent regions of energy have to be distinguished in the one- 
dimensional Anderson model. For energies 2 — > D^^^ 
close to the band edge, corresponding to relatively strong dis- 
order, the RPA fails and the distribution function deviates from 
the strict log-normal form |27|. RPA fails also for energies 
\E\ < close to the band-center |34|, and the distribution 
function again deviates from the strict log-normal form 1261 . in 
generaUzation of the Kappus-Wegner anomaly of Zioc at i5 = 



0.25 



0.15 




0.05 



Fig. 4. The probability distribution function P(ln v) from a numerical 
simulation (data points) in the Anderson model with L — Zioc/2, E = 
— 1, ^pcrt = 24 is compared to a normal distribution (dashed line) with 
the same mean and variance as the data, and a generalized normal 
distribution from the Pearson system (solid line), Eq. J6j, where the 
four free parameters are determined from the asymptotic values of the 
first four cumulants. 



I24I35I36I . For other energies inside the band, the RPA is jus- 
tified, and SPS holds, for weak disorder. 

Offsite disorder adds another means to depart from SPS at 
the band center, since the balance parameter S interpolates be- 
tween the Kappus-Wegner anomaly at (5 = 1 {Dt = 0) and the 
Dyson singularity at 5 = — 1 (Dy = 0), which results in to- 
tal delocalization 1 37 38 '39 40 41 1. In the vicinity of the band 
center, we derive the following Fokker-Planck equation for the 
joint distribution function P{u, a; i) of u ~ — In 5 ~ — In ly, 
and the phase a for reflection from the system: 



da.i^S2o.-e) + dlil + Scl) 



P, (7) 



where = s'mx, = cosx, e — E/{Dv + Dt), and d 
denotes partial derivatives. For 5 = \, this equation has been 
used to study the wavefunction statistics at the Kappus-Wegner 
anomaly |27|. For 5 = —1, E = {), one recovers the dereal- 
ization at the Dyson anomaly. At the balance point of onsite 
and offsite disorder i5 = 0, the variables u and a decouple. For 
L ~> 00 the reflection phase a becomes completely random, 
while P{u) precisely takes the form of SPS, with asymptotic 
slopes C2 — 2ci and c„ = for n > 3. Hence, somewhat 
surprisingly, we find that RPA and SPS hold true by a partic- 
ularly simple mechanism just in between the two abovemen- 
tioned anomalies. 

Away from the novel SPS point (5 = 0, but for disorder 
still small, the asymptotic behavior of the distribution function 
can be analyzed by introducing into Eq. the large-deviation 
ansatz 

/ + ioo jy- 00 
-i.^exp[c(C)i-e«]/(a;0. (8) 
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Fig. 6. Dependence of the parameters Xg""' and X^""' on Xj""' for various ways to depart from SPS [SPS conditions are indicated by the dot 
at coordinates (2,0)]. (a) The balance of weak onsite and offsite disorder is changed at the center of the band, E — 0. The curve is the prediction 
of the analytical theory, while data points are obtained by numerical simulations, (b) Analytical results as energy is varied at a fixed strength of 
onsite disorder, with /pert 3> 1 but /port/Ap arbitrary, around the band center (E — 0) and around the band edge (\E\ ~ 2). (c) Numerics and 
linear fits for variable strength of onsite disorder at three different values of the energy, (d) Comparison of the curves in (a-c). 



Here c(^) = X^n C"'^"/'^' is '^he generating function of the 
slopes of the cumulants, see Eqs. ([0 and (|2jl, and /(a; £,) has 
to be periodic and normaUzable with respect to a. We arrive at 
a differential equation 



do.{-s2o. - e) + dl{l + del) 



Dv + Dt 

+ -(1 + 6c2a - 2daSs2a) + y (1 ^ <5c2q) 



(9) 



in which the slope-generating function c(^) appears as an eigen- 
value, while appears as an eigenfunction. The slopes 
c„ now can be calculated iteratively by expanding c(^) and 
/(a; ^) order by orderin ^, following Refs. |26 27|. Awayfrom 
the SPS point J = but for \E\ < Dy + Dt, the slopes c„ take 
finite values, in compliance with Eq. Q. Our analytical results 
are confirmed by numerical computations in Fig.|5l 



3.2.2 Constant offsets d„ 



The offsets dl"' « CT^ [L = 0) can be calculated by express- 
ing the local density of states v{L) in terms of the reflection 
coefficients tr (ri) from the segment of the wire to the right 
(left) of the point L at which v is calculated I30l . 



v{L) = Re 



(l + rL)(l + rR) 



(10) 



where we normalized {v{L)) 



1 (which amounts to multi- 
and the 



plication by a constant factor 7r\/4 — E'^). For L 
opening of the wire oriented to the left, tl = because there 
is no reflection from the opening, and = exp(za), where a 



is the phase of reflection from the semi-infinite system. Hence, 
the numbers 



« C(-)(L = 0) = (([-ln(l +cosa)n) 



(11) 



characterize the distribution of the reflection phase a of the 
semi-infinite system 1 30 1, and allow to assess the validity of the 



RPA, which predicts d^^^ = In 2, d^^^ = 7rV3, d^^^ = 12C(3) 
(with the Riemann C, function), d^^^ — 147r'*/15. 

cIi\l = 0) vanish independently of 



The offsets d^n^ 



the RR4 since in terms of the reflection matrices introduced 
above 

1 + ri^r-R, 



v{L) = Rc 



1 - TLrR 



(12) 



and hence !>(0) = 1. 

The offsets dn'' are obtained by considering the composi- 
tion law tR+L = iR(l — t'lJ'rJ^^^l for the series transmis- 
sion through two long segments R and L. The reflection coef- 
ficients now are equivalent phase factors tr l = cxp(iQ;R,L)- 
We equate the cumulants of both sides and insert the asymp- 

toticsO. The constant offsets follow as dj'i'' = (-l)"(({ln[2- 
2cos(Q!R+aL)]}"))-IntheRPA,4^^ = Oanddjf^ = {-lYd^n^ 
This is clearly displayed in Fig. ^ Beyond the RPA, the dn^ 
and dn^ contain equivalent information on the reflection-phase 
distribution function P{a), but no longer are simply related. 

3.3 Independence of the parameters 

Now we turn to the question of the mutual independence of the 
parameters c„, as we violate the conditions for SPS. 

A convenient set of parameters beyond the SPS quantity Ci 
is formed by the ratios X„ = C„/Ci, which asymptotically 



6 



H. Schomerus, M. Titov: Short-distance wavefunction statistics in one-dimensional Anderson localization 




-1 1 2 4 6 8 10 



8 \E\l{Dy+D,) 

Fig. 5. Cumulants Cn (in units of "^'^ — 2L/Zpcrt) in the asymptotic 
limit L 2> ^loc, as a function of the balance parameter 5 for energy 
E = Q (left panels), and as a function of E for 5 = and 5 = 1 
(right panels). Results of numerical simulations with /pert = 300 are 
compared to the predictions of the analytical theory. 



acquire the constant values 

lim X„ = - c„/ci. (13) 

L— !-oo 

In RPA and SPS, only one effective parameter Ci survives 
since = = 2ci^^ = 2cl''\ i.e., xi°°^ = 2, and 

moreover Cn'^ — cif* — for n > 3, which gives a picture 
consistent with SPS. However, beyond this approximation the 
cumulants C„ with n > 3 generally may increase linearly with 
L, and hence can be of the same order as C\ and C2, such that 
all are of order unity. Notice that the asymptotic value 

xi""-* is well approximated by Xn^ even for L < Zioc, since 
the cumulants Cn^ are linear already for small L and the off- 
sets d'n^ vanish. 

In Fig.|6lwe plot the asymptotic ratios of cumulants x'^'' 
and X^^^ as function of X^^^^ , while we vary: 

(a) the balance parameter (5 at i? = (a i) from to 1 and 
(a ii) from to —1; 

(b) energy for fixed onsite disorder (b i) around = and 
(b ii) around \E\ — 2; and 

(c) the disorder strength from Zpcrt — 300 to /pert = 1-5 at 
(5 = 1 for the three values of energy (c i) i? = 0, (c ii) E = —1, 
and (c iii) £: = -1.95. 



In the cases (a) and (b) we show the results of the analyti- 
cal procedure described above, while for (c) we show the result 
of the numerical simulations. For illustration of the predictive 
power of the theory presented in the previous Sec. 13.21 numer- 
ical results are also displayed for case (a). 

Of particular interest is the curve for case (b ii), for energies 
close to the band edge, which also applies to strong disorder, 
JJ2/3 > 2 - |27|. [See also the data points for case (c iii).] 

In this case the curves Xn°^^ {X2°°^) depart from the seemingly 
unique functional behavior obtained in the other cases. Hence, 
we are led to conclude that at least for sufficiently strong dis- 
order X^°°^ and x'l^^^ are not uniquely determined by X2°°\ 
Since the SPS quantity Ci is always an independent scaling 
parameter, altogether more than two quantities are needed to 
characterize the distribution function P{i^) for short distances 
L < ^loc (where the central-limit theorem, and hence TPS, not 
yet applies). 



4 Conclusions 

We observed that the short-distance statistics of localized wave- 
functions inside a long one-dimensional disordered system can 
be recovered from the long-distance statistics, but in general 
are characterized by more than the two parameters (a mean Ci 
and a variance C2) that suffice to describe the long-distance 
statistics themselves. These additional parameters have been 
obtained from the higher cumulants C„ of lii i', where ly is the 
local density in a semi-infinite system. The additional param- 
eters in the C„ can be neglected when considering the case 
of weak disorder and generic energies within the band: Then 
Ci ~ C2/2 (X L, and also C„ = 0(L°) for n > 3 take uni- 
versal values, which results in a picture consistent with single- 
parameter scaling even in the short-distance wavefunction statis- 
tics. 

With three-dimensional systems in mind, it would be desir- 
able to investigate the relation of the parameters from large de- 
viation statistics to the scaling parameters at the metal-insulator 
transition, which may be established by multi-fractal analysis 
when this transition is approached from the localized regime. 

Another potential application of the higher cumulants is to 
use them for detecting spatial correlations in the disorder, since 
the higher cumulants are sensitive to higher-point wavefunction 
correlations. This offers a natural extension of a previous in- 
vestigation 1 42 1, which demonstrated that deviations from ran- 
domness due to spatial three-point correlations (such as dis- 
played by a folded Fibonacci sequence) cannot be detected by 
the conventional wavefunction statistics. 
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